Multiqubit symmetric states with high geometric entanglement 
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We propose a detailed study of the geometric entanglement properties of pure symmetric A'^-qubit 
states, focusing more particularly on the identification of symmetric states with a high geometric 
entanglement and how their entanglement behaves asymptotically for large A'^. We show that much 
higher geometric entanglement with improved asymptotical behavior can be obtained in compari- 
son with the highly entangled balanced Dicke states studied previously. We also derive an upper 
bound for the geometric measure of entanglement of symmetric states. The connection with the 
quantumness of a state is discussed. 
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I. INTRODUCTION 

Entanglement of compound systems is known to be a 
key resource in quantum information allowing non-local 
correlations T|. In many practical applications it is of 
fundamental importance to know whether a state is en- 
tangled or not. Very often, however, this information is 
not sufficient and it is also required to know how much a 
state is entangled. A useful tool to quantify the amount 
of entanglement of a state is given by the so-called en- 
tanglement measures A widely used such measure is 
provided by the geometric measure of entanglement that 
is defined for a pure state \ip) by [1] 



Eaim = 1 
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where the maximum is taken over all separable pure 
states. The extension to mixed states is done via the 
usual convex- roof construction |2| . The geometric mea- 
sure fulfills all the requirements of an entanglement 
monotone 31 and is comprised between and 1. While 
the lower value is realized for separable states, the upper 
value is never reached because that would imply having 
a state orthogonal to all separable states and this is sim- 
ply impossible. The geometric measure of entanglement 
has been shown to have several operational interpreta- 
tions, on one side in terms of the utility of a state as 
an initial state for the Grover algorithm Q, and on the 
other side in terms of the difficulty of state discrimina- 
tion under local operations and classical communications 
(LOCC) [a. It is also related to optimal entanglement 
witnesses [3, @ and has nice applications in many-body 
physics and condensed matter systems [7l-[To|. 

Despite its usefulness, the explicit value of the geomet- 
ric measure of entanglement has only been derived, so far, 
for a limited number of entangled states, such as A^-qubit 
Greenberger-Horne-Zeilinger (GHZ) states Q, Dicke 
states 0], generalized W^-states graph states [l^ . 



or other typical states with given symmetry proper- 
ties [l3l - [l6| . The geometric measure remains unknown 
for most of the multipartite states simply because of 
the definition that involves an optimization procedure 
over the class of separable states and this represents a 
formidable task in the general case even with numerical 
approaches. A very recent study showed, however, that 
this task is drastically simplified in the case of states that 
are symmetric under any permutation of the parties fv^ . 
In that case the optimization can be done on the re- 
stricted set of symmetric separable states |0, . . . , 0). This 
represents a huge simplification in the calculation of the 
geometric measure since the number of parameters in- 
volved in the optimization procedure does not depend on 
the number of parties anymore. For qubit systems, only 
two real parameters are required. In this paper we exploit 
this great simplification to propose a detailed study of the 
geometric entanglement properties of symmetric A'^-qubit 
states. We focus more particularly on the identification 
of symmetric states with a high geometric entanglement 
measure and how their entanglement behaves asymptot- 
ically for large A^. We show that much higher geometric 
entanglement with improved asymptotical behavior can 
be obtained in comparison with the highly entangled bal- 
anced Dicke states We also derive an upper bound 
for the geometric measure of entanglement of symmetric 
states, that is significantly tighter than the one known to 
hold for any state [l^ . We finally make a connection be- 
tween the most geometrically entangled states and states 
having the highest quantumness as defined in Ref. [l9i |. 

The paper is organized as follows. In Sec. |lT] the use- 
ful representations of A^-qubit symmetric states are given 
to get a simplified expression of the geometric measure 
of entanglement for that class of states. In Sec. IIIII we 
derive the announced upper bound. In Sec. IIVI we in- 
vestigate various configurations of highly geometrically 
entangled symmetric states. Their asymptotic behavior 



2 



with respect to the number N of qubits is analyzed. The 
connection of the geometric entanglement with the quan- 
tumness of a state is given in Sec. |Vl We finally draw 
conclusion in Sec. I VII 



II. 



GEOMETRIC ENTANGLEMENT FOR 
SYMMETRIC STATES 



K being the polynomial degree and the remaining ai 
equal to 1 [13]. 

The Majorana representation is particularly useful to 
get a simplified expression of the geometric measure of 
entanglement for any symmetric iV-qubit state {ips}- In- 
serting Eq. ^ into Eq. ^ and considering that the max- 
imization is only required over the set of symmetric sep- 
arable states \4>,. . . yields immediately 



In an A^-qubit system, an arbitrary symmetric state 
is conveniently represented in either the Majorana [l^ 
or the Dicke state [2l| representation. In the Majorana 
representation, any symmetric state {ips) is univocally 
specified by N single-qubit states = ai\0) + /3i\l) 

= 1, . . . , N) through a sum over all permutations a of 
the qubits : 



N 



■ , 0cr(Af)) 



(2) 



where 



, ba-(i)i ■ ■ ■ J 4'a-{N)) denotes the product state 
\4>a(i)} (X" • • . (8) |0cr(Ar)) and A/" is a normalization pref- 
actor. In the Bloch sphere picture [23], the single-qubit 
states \(f>i) can be represented as points on the unit sphere 
labelled by two angles {9i,(pi) with Ui — cos(6'i/2) and 
j3i = e*"^' sin {6i/2). Any arrangement of N points on the 
unit sphere thus defines univocally a symmetric state in 
the Majorana representation. 

A particularly simple yet important example is given 
by the symmetric Dicke states with k excitations 



\DNik)) = 



1 



:^|0_31^) 



(3) 



N-k 



with the binomial coefficient of N and k. In the 
Majorana representation, the Dicke states correspond to 
N — k points at the North pole {6 — 0) and k points at 
the South pole of the Bloch sphere. 

In the Dicke state representation, any symmetric state 
\ips) is merely expanded in the orthonormal basis formed 
by the N + 1 Dicke states in the symmetric subspace : 



N 



IV's)=AA^4|i?iv(fc)), 



(4) 



k=0 



with dk {k = 0, . . . , A'^) the complex expansion coeffi- 
cients. 

The Dicke state representation of a symmetric state 
written in the form of Eq. ^ is obtained through the 
relation [23|] 



di. 
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V(JV)- 



(5) 



Inversely, the Majorana representation of a symmetric 
state expressed in the Dicke state basis is obtained with 
N single-qubit \(j)i) states defined by ai/ Pi equal to the K 
roots of the polynomial P{z) = Y.k i-^fiCNf^dkZ^ ^ 



i?G(l^s)) = l-AA27V!2maxTT|((/., 

\4>) 



(6) 



The geometric measure of entanglement of the A^-qubit 
GHZ states IGHZat) = (|0...0) -f |1...1))/V2 is equal 
to 1 /2 regardless the number of qubits [3] . For the sym- 
metric Dicke states, it reads [1] 



EG{\DM{k))) = l-C_ 



N 



N ~k 
N 



N-k 



(J) 



It is maximal for a balanced number of excitations, i.e., 
ioi k = kN with kN closest to N/2. In that case we get 
the asymptotic behavior for large N 



(8) 



If the right-hand side of Eq. (jH)) converges to 1 with the 
number of qubits, this convergence can be qualified as 
quite slow since the difference with one only decreases as 
the inverse of the square root of N. In the next section, 
we derive an upper bound of the geometric measure of 
entanglement that allows one to expect symmetric states 
with a geometric entanglement converging much faster to 
one with the number of qubits. Various configurations of 
symmetric states are identified with that behavior. 



III. GEOMETRIC ENTANGLEMENT UPPER 
BOUND 

The geometric measure of entanglement of an A^-qubit 
symmetric state {ips) necessarily verifies 



EgH^s)) < 1 



1 



A^+ 1 



(9) 



Equation Q is immediately obtained from the rep- 
resentation of the identity in the (A^ + l)-dimensional 
subspace of symmetric states as a combination of all 
projectors onto the symmetric separable states |$) = 



l-N+l 



N + 1 
47r 



sin 0d0 

"'0 



(10) 



where the single-qubit state 10) is parametrised by the 
two angles {6,Lp) according to |0) = cos (0/2) |0) -t- 
e"^ sin (0/2) |1). Since obviously the maximal overlap of 
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\ips) with a symmetric separable state is always larger 
than the average overlap over all symmetric separable 
states 1$), we have 



max|(?As|*)|' > ^ 
|4>) 47r 



sm0d0 / dip\{i;s\^)\^ 
Jo 



(11) 
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> 



N+l' 



which proves the upper bound (|9l). It is interesting to 
note that this upper bound is significantly tighter than 
the one known to hold for any state \tjj) beyond the sym- 
metric subspace given by [1^ 



Earn) < 1 



1 

2^ 



(12) 
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(.Pi, <^ij 
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(0,0), (^0,0), (eo,27r/3), (eo,47r/3) with 
6^0 ~ arccos(— 1/3) 


2/3=0.6666. . . 


5 


(0,0), (^0,0), (00, V2), (00, tt), 
(6io,37r/2) with do = 1.8737... 


0.7011... 


6 


(0,0), (7r,0), (V2,0), (V2,V2), 
(V2,^), (^/2,3V2) 


7/9=0.7777. . . 



TABLE I: Most entangled symmetric states defined by their 
set of angles {6i,ipi). For N = 5, cos0o = a; is given by the 
largest real root of the polynomial 101-|-362a;-|-308a;^+894a;^ + 
670a;'' + 894a;^ + 308a;*' + 362a; + 101a;* , and the corresponding 
geometric entanglement Eg = y is given by the smallest real 
root of the polynomial 25-475y + 7630y^ - 18980y^ + 12824y*. 





IV. HIGH GEOMETRIC ENTANGLEMENT 
CONFIGURATIONS 

A. Highest geometric entanglement configurations 

Although the calculation of the geometric measure of 
entanglement is greatly simplified for the particular case 
of A^-qubit symmetric states, the quest of the states in 
the symmetric subspace having the highest geometric en- 
tanglement remains a task that cannot be solved ana- 
lytically in the general N case. Even numerically the 
task gets very quickly extremely involved for increasing 
N values. Here we report results up to iV = 6 [26]. For 
N — 2, the geometric measure of entanglement of a state 
is equal to the minimal Schmidt coefficient 0] and the 
most entangled symmetric state is thus given by the Bell 
state (|00) -I- |ll))/y2 with Eg = 1/2. For N = 3Mt 
is given by the W state 1153(1)) with Eg = 5/9 
For iV = 4 to 6, the Bloch sphere Majorana represen- 
tation points of the most entangled symmetric states we 
identified numerically are reported in Table U and further 
illustrated in Fig [TJ It can be noticed from the figure 
that these states are characterized by N distinct points 
on the Bloch sphere with a large spread, similarly as N 
equal electrical charges tend to place as far as possible 
from each other when they are constrained to a con- 
ducting sphere (Thomson problem). For A'^ = 4 and 6, 
the points are vertices of two platonic regular polyhedra, 
namely the tetrahedron and the octahedron, respectively. 
For TV = 5, the point configuration is a square pyramid. 
The Dicke state representations of these most entangled 
symmetric states read \/l/3(|-D4(0)) -I- -y/^ |Z?4(3))) for 
iV = 4, ^1-^2 1^)5(0)) ~^\D5{^)) with ^ = 0.8373... 
for = 5, and ./lj2 {\De{l)) + \Dg{5))) for iV = 6. 



FIG. 1: Majorana representation of the maximally entangled 
symmetric states for A'^ = 4 — 6 (from left to right). The blue 
dots are the Majorana representation points of the maximally 
entangled symmetric states. The red dots correspond to all 
closest separable states. There are as many of them as the 
number of polyhedron faces. 



B. Coulomb and related arrangements 

Although similar, the Thomson problem mentioned 
above remains distinct from the quest of maximal entan- 
glement configurations since it requires to find A'^ charge 
positions on a sphere with the different constraint of 
minimizing the total electrostatic energy 

N N 

^ = EE^- (13) 

These configurations of points, which we call the 
Coulomb arrangements, are given in Ref. p7j for A^ up 
to 130. For A^ = 4 to 6, they coincide with the Majorana 
representation points of the most geometrically entangled 
symmetric states, except for N = 5 where the Coulomb 
arrangement is a triangular dipyramid corresponding to 
a geometric entanglement of 0.6875. The Coulomb ar- 
rangement is thus proven not to provide the best point 
configuration for the highest geometric entanglement of 
symmetric states. It is nevertheless expected to provide 
high entanglement configurations in view of the large 
spread it leads to. We illustrate in Fig. [5] the geometric 
entanglement of the symmetric states with Majorana rep- 
resentation points distributed according to the Coulomb 
arrangements. These states are denoted by |'0Coui(Af)). 
We find numerically that the entanglement of these states 
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FIG. 2: Geometric entanglement of symmetric states for the 
Coulomb arrangement for A*" up to 110 (blue circles). Equa- 
tion (|14p is displayed in red. Green triangles correspond to 
the entanglement of symmetric Dicke states with a balanced 
number of excitations given by Eq. The grey shaded area 
shows the domain ruled out by the upper bound 
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FIG. 3: Geometric entanglement of symmetric states for dif- 
ferent arrangements : Coulomb (blue circles), Tammes (or- 
ange triangles) and covering (red squares). The entanglement 
of Dicke states with a balanced number of excitations is dis- 
played in green (reversed triangles). The grey shaded area 
shows the domain ruled out by the upper bound ((9]). 



scales up to small fluctuations like 



ScdV-CouKA^))) ^ 1 



C 



iV+ 1' 



(14) 



with C ~ 1.71 a numerical constant. It is very interest- 
ing to note that the Coulomb arrangements, though they 
do not necessarily provide the most geometrically en- 
tangled symmetric states, define nevertheless states with 
a very high geometric entanglement close to the upper 
bound ^ and behaving asymptotically with the number 
of qubits in a similar fashion. They can therefore de- 
fine a good alternative strategy to get highly entangled 
symmetric states. 

We have investigated other configurations of Majorana 
representation points by considering other ways of dis- 
tributing as evenly as possible N points on a unit sphere. 
The Tammes and covering problems define such alterna- 
tive ways. In the Tammes problem, one seeks to max- 
imize the minimal distance min^^j jr^ — r^l among the 
N points on the sphere. In the covering problem, one 
rather seeks to minimize maxrgs miuj |r — r^j that repre- 
sents the greatest distance between a point r on the sur- 
face of a unit sphere S and the nearest of the A'^ points. 
Tammes and covering arrangements are also given in 
Ref. j27.] for TV up to 130. For TV = 2 - 6 and iV = 12, 
Coulomb, Tammes and covering arrangements are iden- 
tical |28'] . Figure [3] shows the entanglement relative to 
these different arrangements together with the entangle- 
ment of symmetric Dicke states with a balanced number 
of excitations. Coulomb and covering arrangements give 
very similar results. In contrast, Tammes arrangement 
shows large fiuctuations. The least fiuctuating curve cor- 
responds to Coulomb. 



C. Arrangements with equally weighted 
superpositions of Dicke states 

Equally weighted superpositions of Dicke states may 
exhibit interesting entanglement properties with respect 
to the geometric measure if their relative phases are prop- 
erly chosen. We first consider the states 



N 



k=0 



VN + T 



\DN{k)), 



(15) 



with 7 a real parameter. The geometric entanglement of 
these states as a function of N is displayed in Fig. |4] for 
two values of 7 (2/3 and 1). It is compared to the geo- 
metric entanglement of the Coulomb arrangement. For 
the two values of 7 under consideration, the states (fTSj) 
are clearly less entangled than Coulomb. However, in 
contrast to the latter, the curves show interestingly al- 
most no fluctuations and the entanglement behaves in a 
similar fashion for large N : at least, for N up to 100, 
we have 



1 - 



N + 1' 



(16) 



with D-y a constant depending on 7 (-D2/3 ~ 2.22, Di « 
2.81). 

We obtained numerically that the Majorana represen- 
tation points of the state (jl5|) spread in a kind of spiral 
on the Bloch sphere, with angles {6k, fk) given by 



arccos 1 



N -1 

ipk fa 7(1 — 2k) mod 27r. 



(17) 
(18) 



The corresponding Majorana representation is illustrated 
in Fig. [S] for iV = 400 and 7 = 2/3. The large spread 
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FIG. 4: Geometric entanglement for different sequences of 
symmetric states as a function of the number A'^ of qubits. 
Blue circle : i;G(|V'Coui(A^)>); red triangles : EGi\'4>iii{N)))\ 
green squares : iJodV'iC-^))); orange reverse triangles : 
-^g(I V'2/ 3(-^))) fitted in excellent agreement (black curve) by 
Eq. (1221) . The grey shaded area shows the domain ruled out 
by the upper bound @. 
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FIG. 5: Majorana representation of the state (|15p for 7 = 2/3 
and iV = 400. 



of the points explains the high geometric entanglement 
value obtained for that state. 

We next consider another equally weighted superpo- 
sitions of Dicke states with linearly increasing phases, 
namely 



N 

lin(7V)) = 



k=0 



VN + T 



\DN{k)). 



(19) 



These states have the advantage that their geometric 
entanglement can be analytically derived for large N. 
Their overlap with the symmetric separable state |$) = 
10, . . . , 0) reads explicitly 



The remaining sum is well approximated by replacing the 
binomial distribution with a Gaussian distribution, and 
the sum by an integral. The maximum of the Gaussian 
distribution is situated at A^(l — p), and its variance is 
Np{l ~ p), so as for sufficiently large N the integration 
range can be extended to ±00. We then get 



^87Tpil-p)N 
N +1 



(21) 



where p = cos^ (^/2). The overlap (I2ip is maximum for 
6 = tt/2 and yields the geometric entanglement 



N + 1 



(22) 



Figure m shows that Eq. ([22]) reproduces the exact nu- 
merical data very well, as soon as N is of the order 10. 
Comparison with Eq. ([S]) shows that superposing Dicke 
states with linearly increasing phases reduces the geomet- 
rical entanglement compared to a single Dicke state with 
k equal (or close to) N/2. The fact that a 1/Vn behav- 
ior is obtained for linear phases, but a 1/iV behavior for 
quadratically increasing phases, demonstrates that the 
phases of the Dicke states play a key role in the amount 
of geometric entanglement that can be achieved. 



V. GEOMETRIC ENTANGLEMENT AND 
QUANTUMNESS 

In [l9| a quantumness measure was introduced in or- 
der to quantify how quantum an arbitrary mixed spin-j 
quantum state p is. The measure was defined as the 
Hilbert-Schmidt distance from the state p to the set of 
classical states (which is the convex hull of spin coher- 
ent states). A spin-j state can be seen as a symmetrized 
state of 2j spins i, and thus it can be written under the 
form Eq. ([2]) with N — 2j. A spin-j coherent state is 
the tensor product of 2j identical spin-i coherent states. 
Since any spin-^ pure state is a coherent state, spin-j co- 
herent states are by definition the symmetric separable 
pure states. 

Let p be a generic spin-j density matrix. We define 
the B-quantumness Qb{p) as the distance from p to the 
set C in the vector space of {N -I- 1) x {N + 1) matrices 
equipped with the Bures metric. Namely 



Qb{p) = min Db{p,Pc), 

Pc^C 

where Db is the Bures distance [1^ defined by 



(23) 



|«"(7V)|$)|2^ 



1 



7V + 1 



N 



Y.^C^p^-^^il-pfe^^'^^-'^^^ 

(20) 



k=0 



with (0, Lp) the two angles of the single qubit state 
The maximization over ip yields straightforwardly Lp - 



7- 



Db{pa,Pb) 



2tr 



PAPBy/PA 



(24) 



In other words, we consider a problem analogous to the 
one in [l^ but replace Hilbert-Schmidt distance by Bures 
distance. If p = \4's){4's\ is a pure state density matrix. 
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with l^/jg) a generic spin-j state defined as in Eq. ([2|), the 
expression for the Bures distance simplifies drastically to 

DB{\i^s){'^slpc) - v^2-2v/(V's|pclV's)- (25) 

Any element p G C can be decomposed as a mixture of a 
certain number n of coherent states i.e. 

n 

PC = ^ A, I I (26) 

4=1 

with the Xi positive weights summing up to 1. The scalar 
product appearing in Eq. (j25p then reads 

TL 

{^Ps\pc\iJs)=Y.^^\(^^\'^s)\'- (27) 

4=1 

B-quantumness is obtained by minimizing the quantity 
(1251) over C, which is equivalent to maximizing the scalar 
product in (|27p. The maximum with respect to varia- 
tion of the non-negative Xi occurs when all but one are 
equal to zero. The problem thus reduces to finding the 
largest overlap of the fixed \ijjs) with a coherent state, 
and Eq. becomes 



QB{Hs){i's\) - ./2- 2max|($|V'5)|. (28) 

Thus B-quantumness and geometric entanglement ^ are 
essentially the same quantity. The maximally entan- 
gled symmetric states obtained in Sec. HV] are thus the 
pure states with largest B-quantumness, which we call 
Queens of B-Quantumness (QBQ). As was mentioned 
in p^j for Hilbert-Schmidt quantumness, one can easily 
deduce from the fact that Qb is a convex function that B- 
quantumness reaches its maximum for pure states, thus 
QBQ states can be looked for among pure states. It is 
interesting to compare these QBQ to the most quantum 
states (QQ states) obtained in J^]. For the lowest values 
of N, the QBQ states differ from the QQ states only for 
iV = 3 and 5. 



In the case of two qubits, entropy of entanglement is 
known to be the unique entanglement measure for pure 
states [lO]. Therefore QBQ state for N = 2 coincides 
with the Dicke state \D2{1)) which is also the QQ state 
(see Sec. IIVI and Ref. ^1^]). For more than two qubits 
there is no unique measure of entanglement. Thus one 
should not expect to find a unique measure of quantum- 
ness, and the discrepancy between QBQ and QQ states 
should come as no surprise. 



VI. CONCLUSION 

In summary, we have realized a detailed study of the 
geometrical entanglement properties of symmetric A^- 
qubit states. We have focused on the identification of 
symmetric states with a high geometric entanglement 
and how their entanglement behaves asymptotically for 
large A^. We have shown that much higher geometric en- 
tanglement with improved asymptotical behavior can be 
obtained in comparison with the highly entangled bal- 
anced Dicke states 3]. We have also derived an upper 
bound Eg = 1 — 1/(A^ +1) for the geometric entangle- 
ment in the restricted set of the symmetric states. This 
value is significantly tighter than the one holding for any 
state. The connection with the quantumness of the states 
has been established. 
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